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We analyze the B — > A'^(— > Kix)l + l~ (with / = e,/j,,r) decay in the standard model and two 
new physics scenarios: vector-like quark model and family non-universal Z' model. We derive 
I ' its differential angular distributions, using the recently calculated form factors in the perturbative 

^-H ■ QCD approach. Branching ratios, polarizations, forward-backward asymmetries and transversity 

amplitudes are predicted, from which we find a promising prospective to observe this channel in the 

C^l , future experiment. We update the constraints on effective Wilson coefficients and/or free parameters 

in these two new physics scenarios by making use of the B — > K*l + l~ and b — > sl + l~ experimental 
data. Their impact on B —¥ K^l^l" is subsequently explored and in particular the zero-crossing 
point for the forward-backward asymmetry in these new physics scenarios can sizably deviate from 
the standard model. In addition we also generalize the analysis to a similar mode B s — » /2(1525)(— > 

<N ! K+K-)i+r. 

(—| ! PACS numbers: 13.20.He; 12.39.St 14.40.Be; 

I. INTRODUCTION 

Discoveries of new degrees of freedom at TeV energy scale, with contributions to our understanding of the origin 
£SJ . of the electroweak symmetry breaking, can proceed in two different ways. One is a direct search of the Higgs boson, 
the last piece to complete the standard model (SM), and particles beyond the SM, to establish new physics (NP) 
theories. The other effort, already ongoing, is to investigate processes in which SM is tested with higher experimental 
and theoretical precision. Among the latter category, rare B decays are among ideal probes. Besides constraints on 
£NJ ' the Cabibbo-Kobayashi-Maskawa (CKM) matrix including apexs and angles of the unitary triangle, which have been 
. contributed by semileptonic b — > u/c and nonleptonic B decays respectively, the electroweak interaction structure can 
t— I ' also be probed by, for instance, the b —> s-f and b —> sl + l~ modes which are induced by loop effects in the SM and 
therefore sensitive to the NP interactions. 

Unlike b — > S7 and B — > K*"/ that has only limited physical observables, b — > sl + l~ especially B — > K*l + l~ , 
^ • with a number of observables accessible, provides a wealth of information of weak interactions, ranging from the 
forward-backward asymmetries (FBAs), isospin symmetries, polarizations to a full angular analysis. The last barrier 
;_i ' to access this mode, the low statistic with a branching faction of the order 10 -6 , is being cleared by the B factories 
and the hadron collider [D-Q . The ongoing LHCb experiment can accumulate 6200 events per nominal running year 
of 2fb^ 1 with y/s = 14 TeV [1], which allows to probe the short-distance physics at an unprecedented level. For 
instance the sensitivity to zero-crossing point of FBAs can be reduced to 0.5GcV 2 and might be further improved as 
0.1 GeV 2 after the upgrade [H]. This provides a good sensitivity to discriminate between the SM and different models 
of new physics. There are also a lot of opportunities on the Super B factory 0]. Because of these virtues, theoretical 
research interests in this mode have exploded and the precision is highly improved, see Refs. [7 - 22 1 for an incomplete 
list. 

Toward the direction to elucidate the electroweak interaction, B — > K*l + l~ and its SU(3)-related mode B s — > tfil + l~ 
are not unique. In this work, we shall point out that B — > if|(1430)^ + ^ and the £? s -countcrpart B s — > (1525)Z + Z — 1 , 



which so far have not been investigated in detail [23j-|26j, are also useful in several aspects. Due to the similarities 
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between K* and K%, all experiment techniques for B — > K*l + l~ are adjustable to 5 -> K^l + l~ ■ The main decay 
product of K2 is a pair of charged kaon and pion which are easily detected on the LHCb. Moreover as we will show 
in the following, based on either a direct computation in the perturbative QCD approach 27| or the implication 
of experimental data on B — > process, the branching ratio (BR) of B — > K^l + l~ is found sizable. Therefore 
thousands of signal events can be accumulated on the LHCb per nominal running year. 

As a consequence of the unitarity of quark mixing matrix, tree level flavor-changing neutral-current (FCNC) is 
forbidden in the SM. When higher order corrections are taken into account, b — > sl + l~ arises from photonic penguin, 
Z penguin and W-box diagram. The large mass scale of virtual states leads to tiny Wilson coefficients in b quark 
decays and thus b sl + l~ would be sensitive to the potential NP effects. In certain NP scenarios, new effective 
operators out of the SM scope can emerge, but in a class of other scenarios, only Wilson coefficients for effective 



operators are modified. Among the latter category, vector-like quark model (VQM) |28l436j | and family non-universal 
Z' model 37-42J are simplest and therefore of theoretical interest. In this work we shall also elaborate the impacts 
of these models on B — s- KZl + l~ . 



The rest of the paper is organized as follows. In Sec. [TIJ we collect the necessary hadronic inputs, namely form 
factors. Sec. II contains the analytic formulas for differential decay distributions and integrated quantities. In Sec. lIVl 
we give a brief overview of two NP models whose effects we will study. Sec. V is our phenomenological analysis: the 
predictions in the SM; update of the constraints on the VQM and Z' model parameters; the NP effect on the physical 
quantities. We conclude in the last section. In the appendix, we give the effective Hamiltonian in the SM and the 
helicity amplitude method. 



II. B -> K 2 FORM FACTORS 

B — > K%l + l~ decay amplitudes contain two separate parts. Short-distance physics, in which contributions at the 
weak scale fiw is calculated by perturbation theory and the evolution between mw and b quark mass scale m& is 
organized by the renormalization group. These degrees of freedom are incorporated into Wilson coefficients and the 
obtained effective Hamiltonian responsible for b — > sl + l~ in the appendix A. The low-energy effect characterizes the 
long-distance physics and will be parameterized by hadronic matrix elements of effective operators, which are usually 
reduced to heavy-to-light form factors in scmilcptonic B decays. 

The spin-2 polarization tensor, which satisfies e^P^ = with Pi being the momentum, is symmetric and traceless. 
It can be constructed via the spin-I polarization vector e: 

e^(±2) = e M (±)e„(±), e^(±l) = -^[e M (±)e„(0) + e„(±) C/1 (0)], 

e^(0) = -L[e M (+)e 1/ (-) + e v (+)e„(-)] + ^|e M (0)e„(0). (1) 

In the case of the tensor meson moving on the z axis, the explicit structures of e in the ordinary coordinate frame are 
chosen as 

e„(0) = —(\p KS \,0,0,E Ki ), e /1 (±) = ^(0,Tl,-i,0), (2) 
rriK* V2 

where Ek* and px* is the energy and the momentum magnitude of K\ in B meson rest frame, respectively. In the 
following calculation, it is convenient to introduce a new polarization vector ex 

e T »{h) = — e^{h)P v B , (3) 

which satisfies 

£7> (±2)=0, e T „(±l) = — -U(0).P B e M (±), e TM (0) = — J|e(0) • P s ^(0). (4) 

m B V2 rriB V 3 

The contraction is evaluated as e(0) • Ps/ms = \pk* \/^nK* and thus we can see that the new vector ex plays a similar 
role to the ordinary polarization vector e, regardless of the dimensionless constants \pk* \ /m>K* or \J^\pK*\/'mK*- 
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The paramctrization of B — > K 2 form factors is analogous to the B —> K* ones 2514271 143 | 

2V(q 2 ^ 



{K* 2 {P 2 ,e)\s^b\B{P B )) 
(KZ(P 2 ,e)\s^ l5 b\B(P B )) 

(KI(P 2 ,e)\sa^q v b\B(P B )) 
(K;(P 2 ,e)\s<T^ l5 q„b\B(P B )) 



m B + rriK. 



fiupa * p p 



2im K *A (q 2 )^Y i q tl + i(m B + m K *)Ai(q 2 ) 



-*A 2 (q 2 )- 



e T • q 



•K 



m B + rriK* 
-2iT x {q 2 )e^ ^ v P Bp P 2a , 

T 2 (q 2 ) (m 2 B - m 2 K ,)e* Tll - e* T ■ qP^ + T 3 (q 2 )e* T ■ q 



'A.? 



(5) 



where q = P B — P 2 , P = P B + P 2 . We also have the relation 2rriK*Ao(0) = (m B + rriK*)Ai(Q) — (m B — rriK*)A 2 (0) 
in order to smear the pole at q 2 = 0. 

Using the newly-studied light-cone distribution amplitudes 44], we have computed B — > K 2 form factors [2t| in 
the perturbative QCD approach (PQCD) (45j . At the leading power, our predictions are found to obey the nontrivial 
relations derived from the large energy symmetry. This consistence may imply that the PQCD results for the form 
factors are reliable and therefore suitable for the study of the semilcptonic B decays. The recent computation in 
light-cone QCD sum rules |43| is also consistent with ours. Results in the light-cone sum rules in conjunction with 
-B-meson wave functions [461 ] . however, are too large and thus not favored by the B — > K 2 "f data. In our work the 
B —> K 2 form factors are g 2 -distributed as 27 1 



F(q 2 ) 



F(0) 



(1 - q 2 /m 2 B )(l - a(q 2 /m 2 B ) + b{q 2 /m 2 B ) 2 ) ' 



(6) 



where F denotes a generic form factor among A$ : A±, V, Ti_3. Neglecting higher power corrections, A 2 is related to 
Aq and Ai by 



A 2 (q 2 ) 



5 f - [(m B +m Ki )A 1 {q 2 ) - 2m Ki A (q 2 )] 

m B q 



(7) 



Numerical results for the B — > K 2 and B s /2(1525) form factors at maximally recoil point and the two fitted 
parameters a, b are collected in table HI The two kinds of errors are from: decay constants of B meson and shape 
parameter uj b ] Aqcd, the scales ts and the threshold resummation parameter c [271 ] . 



III. DIFFERENTIAL DECAY DISTRIBUTIONS AND SPIN AMPLITUDES 

In this section, we will discuss the kinematics of the quasi four-body decay B — > K 2 {-^ Kir)l + l~ , define angular 
observables and collect the explicit formulas of helicity amplitudes and/or transversity amplitudes. 



A. Differential decay distribution 

At the quark level, the decay amplitude for b — > sl + l~ is expressed as 

\AI1 , 1+1- \ Gf C( em I Cg + Cio r _, , ,j n Cg — Cio r _,, r y ;1 

M{b^sl + l ) = — V tb V t * s x [sb} v _ A [ll\ v+A + [sb] v _ A [ll] v _ A 

\2 7T \ 4 4 

+C 7L m b [sia^(l + 7 5 )&]\ x [l-y"l\ + C 7R m b [sia^{l - 7s)&]^ x , (8) 

where Cjl = Cj and Cm = ^-CVl in the SM. Sandwiching Eq. ([8]) between the initial and final states and replacing 
the spinor product [sb] by hadronic matrix elements, one obtains the decay amplitude for hadronic B process. For 
the process under scrutiny in this work, the decay observed in experiment is actually B — > K 2 (— > Kir)l + l~ which is 
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TABLE I: B — > A'J and B a — » /2(1525) form factors in the PQCD approach. F(0) denotes results at g 2 = point while a, 6 
are the parameters in the parametrization shown in Eq. ©. The two kinds of errors are from: decay constants of B meson 
and shape parameter ut>; Aqcd, factorization scales ts and the threshold resummation parameter c. 
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FIG. 1: Kinematics variables in the B — > A|(— s> K~ ir + )l + l~ process. The moving direction of AJ in B rest frame is chosen 
as the z axis. The polar angle Ok is defined as the angle between the flight direction of K~ (fj,~) and the z axis in the AJ 
(lepton pair) rest frame. The convention also applies to B s — > K + K~)l + l~ transition. 



a quasi four-body decay. The convention on the kinematics is illustrated in Fig. [T] The moving direction of K% in B 
meson rest frame is chosen as z axis. The polar angle 8k (61) is defined as the angle between the flight direction of 
K~ ) and the z axis in K\ (lepton pair) rest frame. cf> is the angle defined by decay planes of AT^ arm the lepton 
pair. 

Using the technique of helicity amplitudes described in the appendix B, we obtain the partial decay width 

d 4 r ^ 

z\Mb\\ (9) 



dq 2 d cos Bud cos 8id<p 8 
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with the mass correction factor /3t = y/T— Amf/q 2 . The function |.M,b| 2 is decomposed into 11 terms 

\M B \ 2 = UlC 2 + 2IIS 2 + (1 C 2 C 2 + 2I S 2 S 2 ) cos(26» ; ) + 2I 3 S 2 sin 2 0, cos(20) + 2^/ 4 C5sin(2^) cos< 
+2v / 2I 5 CS , sin(6'0 cos(/> + 2I 6 S 2 cos6» z + 2\/2I 7 CS sin(^) sm.(j) 
+2\/2I 8 CS sin(2<9 ; ) sin + 2/gS 2 sin 2 6 t sin(2<£) , 

with the angular coefficients 

9 9 
■i ll i ii 

II = (\A LQ \ 2 + \A m \ 2 ) + 8-LRc{A L0 A R0 ] + 4-f \A t \ 2 , 



(10) 



-ff(\A L0 \ 2 + \A R0 \ 2 ), 

lfi(\A L± \ 2 + \A L]] \ 2 + \A RX \ 2 + \A Rll \ 2 ), 
1 



II = 7 [\A L± \ 2 + \A Lll \ 2 + \A R± \ 2 + \A Rll \ 2 ](l 

I c 2 = 

h = 

h = 

h = 

h = 



4m 



3<j 2 



Ami 



Re[A L± A* R± + A L{l A* Rll ], 



P 2 (\A L± \ 2 -\A Lll \ 2 + \A R± \ 2 ~\A m \ 2 ), 



-.pf[Rc(A L0 A* L ,,) + Re(A R0 A*j 



h = \/2pi[Re(A L0 Al ± ) - Re(A R0 A RX )}, 



2j3i{Re(A L]l A* L± ) - Re^,^)], I 7 = V2f3,[lm(A LQ A* Lll ) - Im^o^,)], 
-±=p 2 lIin(A L0 A* L± ) + lm(A R0 A* RX )], 7 9 = /3 2 [lm(A Lll A* L± ) +lTn(A Rll A* R± )}. 



(11) 



C = C(K 2 ) and S = S(K 2 ) for B —> K 2 l + l~ . Without higher order QCD corrections, I 7 is zero and are tiny 

in the SM and the reason is that only Cg has an imaginary part. In this sense these coefficients can be chosen as an 
ideal window to probe new physics signals. 

The amplitudes Ai arc generated from the hadronic B K 2 V amplitudes Hi through Ai = 



N K 



Vx 



2 \/6m B m K * 2m K * yfq 2 



+2m b (C 7L - C 7R )[(m 2 B + 3m 2 - q 2 )T 2 



(Cg - C 10 )[(m 2 B - m 2 K * - q 2 )(m B + m K *)Ai 
A 



m B + m K f 



-A, 



72 Z^ T 3l 



m B ~ m K> 



A 



L± 



N K . 



Vx 



(Cg - C w )[(m B + m K ')Ai =p 



m B + m K * 



V] 



- V8m B mK. 

2m b (C 7L + C 7R ) r- 2m b (C 7L -C 7R ) 2 2 
- (iVATi) + 5 (m B - m K ,)T- 



A 



LI 



r 

n k;^^ (C9 - C W )^=A , 

\J t om B m K * V q 1 

1/2 



(12) 



with Nk* = [ 3 .^il"l° a \VtbV t * s \ 2 q 2 X^ 2 (l - ^) B(K* -> Kir)} 1 / 2 . For convenience, we have introduced transver- 
sity amplitudes as 



.4 



i-L/ll 



A L1 _ = 



A 



-j=(A l+ tA l _), 
-V~2 ^ N KS 



m B m K * 



(Cg 



VXV + 2m b (C 7L + C 7R)v - XTi 



m B + rriK. 



V2^-^ Nkx 



m B m K * 



(Cg - C w )(m B + m K *)Ai + 



2m b (C 7L - C 7R ) 



(m 2 B - m 2 K , )T 2 



(13) 
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and the right-handed decay amplitudes are similar 

A Ri = A L i\c 10 ->-C 10 - 
The combination of the timelike decay amplitude is used in the differential distribution 



A t = A Rt - Aja = 2Nk; 



y/6m 



r 



(14) 
(15) 



B. Dilepton spectrum distribution 



Integrating out the angles 0i , 9k and (f>, we obtain the dilepton mass spectrum 



dT 

dq 2 



1 



(3/ 1 c + 6J 1 s -/ 2 c -2/|) 



and its expression in the massless limit 



\A 



m 



(16) 



(17) 



with i = 0, ±1 or i — 0,1., ||. After some manipulations in the appendix, the correspondence of the above equations 
and Eq. (|20|) with results in Ref. [25j can be shown. 



C. Polarization distribution 



The longitudinal polarization distribution for B — > K 2 l + l is defined as 

df L _ dT Q I dT 37f - I r i 



dq 2 dq 2 I dq 2 3/f + 6Jf - If - 2i| 



(18) 



in which 4% can be reduced into If in the case of m; = since If = — 1%. The integrated polarization fraction is 



given as 



(19) 



D. Forward-backward asymmetry 



The differential forward-backward asymmetry of B — > K 2 l + l is defined by 



dA 



FB 



dq 2 



dcos I 



d 2 T 



= 5/., 



while the normalized differential FBA is given by 



dA 



dq- 

In the massless limit, we have 

dA FB Xq 2 G 2 F al m 



dq 2 

dr 
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' dq 2 d cos 8 \ 4 



3/ B 



3/f + 6q - I r i - 21% 
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FribimB + rriK* ) 
C 9 C W A{V + C W (C 7L + C 7R ) — = *-A Y T x 



+C 10 (C 7L - C 7R ) by * ^1T 2 V 



In the SM where C 7R is small, the zero-crossing point sq of FBAs is determined by the equation 

mh(mR + mK*) mhfmB — tor-*) 
CbAi(a )V(«o) + C 7L -±^ ^Ai(5 )Ti(«o) + C 7L -±-^ ^t 2 ( So )V( S o) = 0. 



•so 



•so 



(20) 



(21) 



(22) 



(23) 
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E. Spin amplitudes and transverse asymmetries 

Using the above helicity/spin amplitudes, it is also possible to construct several useful quantities which arc ratios 
of different amplitudes. The following ones, widely studied in the B — > K* case, are stable against the uncertainties 
from hadronic form factors 



with the notation 



-2Re(i4||Al) 

" r_ + r+ - |Al| 2 + I4iI 2! 

^ r \A i I 2 — I /4m I 2 



(2) = \A^-\A^ 
\A ± \i + \A r ^ 



aP 



II 

Alo^2|| + ^H0^fl|| I 



V\M 2 \a±\ 



2 



4 (4) I^W^Ll -^flo^j_l , 94 x 



= A Ll A* L] + A m A* Ry (25) 



Due to the hierarchy in the SM T_ ^s> T + , A>P is close to 1 and therefore its deviation from 1 is more useful to reflect 
the size of the NP effects. 



IV. TWO NP MODELS 

The b —> sl + l~ has a small branching fraction since the SM is lack of tree level FCNC. It is not necessarily the same 
in extensions. In this section we will briefly give an overview of two NP models, which allow tree-level FCNC. Both of 
these two models, vector-like quark model and family non-universal Z' model, do not introduce new type operators 
but instead modify the Wilson coefficients Cg,CiQ. To achieve this goal, they introduce an SU{2) singlet down-type 
quark or a new gauge boson Z' . 



A. Vector-like quark model: Z-mediated FCNCs 

In the vector-like quark model, the new SU(2)l singlet down quarks Dl and Dr modify the Yukawa interaction 
sector 

C Y = Q L Y D Hd R + h D Q L HD R + m D D L D R + h.c, (26) 

where the flavor indices have been suppressed. Ql (H) is the SU(2) quark (Higgs) doublet, Yd and ho are the 
Yuakwa couplings and mjj is the mass of exotic quark before electroweak symmetry breaking. When the Higgs field 
acquires the vacuum expectation value (VEV), the mass matrix of down type quark becomes 

m d = - - - , (27) 
V | m D J 

which can be diagonalized by two unitary matrices 

mf a = VBmiVg. (28) 



The SM coupling of Z-boson to fcrmions is flavor blind, and the flavor in the process with exchange of Z-boson 
is conserved at tree level. Unlikely although the right-handed sector in the VQM is the same as the SM, the new 
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left-handed quark is SU(2)l singlet, which carries the same hypercharge as right-handed particles. Therefore the 
gauge interactions of left-handed down-type quarks with .Z-boson are given by 

Cz = Ql— g -r-(h -sin 2 Q w Q)$Q L + D L —2— (-sin 2 6 W Q)$D L , (29) 

COS Uw COS Uw 

where g is the coupling constant of SU(2)l, 9\v is the Weinberg's angle, Pr(l) = (1 ±7s)/2. h an d Q are operators 
for the third component of the weak isospin and the electric charge, respectively. 

Since the ratio £d of the coupling constants deviates from unity: £d = — sin 2 OwQd/ — sin 2 OwQf), tree level 
FCNC can be induced after the diagonalization of the down-type quarks. For instance, the interaction for b-s-Z in 
the VQM is given by 

= 9 -^-^s 1 ^P L bZ ll + h.c, (30) 
cos 9w 

where X s b is introduced as the new free parameter: 

A s6 = - l)(V£) sD (V£y bD = |A s6 |exp(i0 s ) . 
Using Eq. (|30l) . the effective Hamiltonian for b — >• sl + l~ mediated by Z-boson is found by 

HLst+i- = ^=-\ sb c s L (- S b) V - A [cl(U) V -A + ci(U) v+A ] ■ (31) 
The Wilson coefficients Cg^io are modified accordingly 

^VLQ _ r SM 47r X sbC S L {c{ + C R ) VLQ _ r SM , ^ A s& c l(4 - C R ) 

ttem Vts V tb a em V ts Vtb 



Making use of the experimental data of b — > sl + l , our previous work 47| has placed a constraint on the new coupling 
constant 

|A s6 | < 1 x 10" 3 , (33) 

but its phase 9 S is less constrained. In the following, we shall see that the constraint can be improved by taking into 
account the experimental data of the exclusive process B — > K*l + l~ . 



B. Family non-universal Z' model 

The SM can be extended by including an additional U{1)' symmetry, and the currents can be given as following in 
a proper gauge basis 

J% = <?' E M4 L Pl + ef R P R ]A, (34) 

i 

where i is the family index and ip labels the fcrmions (up- or down- type quarks, or charged or neutral leptons). 
According to some string construction or GUT models such as Eq, it is possible to have family non- universal Z' 
couplings, namely, even though ef' R are diagonal the gauge couplings are not family universal. After rotating to the 
physical basis, FCNCs generally appear at tree level in both LH and RH sectors. Explicitly, 

= V} L , B^ = V\ R . (35) 

Moreover, these couplings may contain CP-violating phases beyond that of the SM. 
The Lagrangian of Z'bs couplings is given as 

£|cnc = -g'{B^s Lltl b L + B^SR^bR)^ + h.c. (36) 

It contributes to the b — > s£ + £~ decay at tree level with the effective Hamiltonian 

= ^(fasL-rpbL + P* b - S RiMMihin L + pSIr^Zr) , (37) 
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where 

L,R _ g'Mz U L,R / QO \ 

P ff = jM^ B ff (38) 

with the coupling g associated with the SU(2)l group in the SM. In this paper we shall not take the renormalization 
group running effects due to these new contributions into consideration because they are expected to be small. For 
the couplings are all unknown, one can see from Eq. (|37p that there are many free parameters here. For the purpose 
of illustration and to avoid too many free parameters, we put the constraint that the FCNC couplings of the Z' and 
quarks only occur in the left-handed sector. Therefore, p^ b = 0, and the effects of the Z' FCNC currents simply 
modify the Wilson coefficients Cg and Cio in Eq. (|A1|) . We denote these two modified Wilson coefficients by Cg and 
Ci , respectively. More explicitly, 



4t P L M+Pu) r z> _ r , 4tt P L sb {pfi- Pu) 



°9 - c 9 _ T7~T7l ' °io - °io + TTTri ■ 



Compared with the Wilson coefficients in the vector- like quark model in Eq. (f52"j) , we can see that the Z' contribu- 
tions in Eq. (|39p have similar forms and the correspondence lies in the coupling constants 

\ sb c s L ^ P L sin c l L:R -> p^. (40) 

However the number of free parameters is increased from 2 to 4 since c l L R in the VQM is the same as the SM. 

V. PHENOMENOLOGICAL ANALYSIS 

In this section, we will present our theoretical results in the SM, give an update of the constraints in the above two 
NP models and investigate their effects on B — s- K 2 p + p~ and B s —> f 2 p + p~ . For convenience, branching ratios of 
K 2 and f 2 decays into Kix and KK will not be taken into account in the numerical analysis. 

A. SM predictions 



With the B — > K 2 form factors computed in the PQCD approach 27(, the BR, zero-crossing point of FBAs and 
polarization fractions are predicted as 

B(B -> K* 2 p+p-) = (2.5±\*) x 10~ 7 , 

f L (B^K* 2 p+p-) = (66.6 ±0.4)%, 

s (B -> Kln + n~) = (3.49±0.04)GcV 2 , 

B(B^K* 2 t+t-) = (9.6^j) x 10- 10 , 

f L (B -> K*t+t-) = (57.2 ±0.7)%. (41) 

The errors are from the form factors, namely, from the B meson wave functions and the PQCD systematic parameters. 
Most of the uncertainties from form factors will cancel in the polarization fractions and the zero-crossing point so. 
Similarly results for B s — > f 2 l + l~ are given as 



B(B S - 


» f' 2 P + V 


■) 


= (i.8±5;|) x 10- 7 , 




» f 2 P + ^ 


■) 


= (63.2 ±0.7)%, 


s (B s - 


» f 2 P + V 


■) 


= (3.53±0.03)GcV 2 


B(B S - 


+ f- 2 r + r- 


■) 


= (5.8±|;I) x 10- 10 , 


f L (B s - 


+ f 2 r + r- 


) 


= (53.9 ±0.4)%. 



(42) 

We also show the q 2 -dependence of their differential branching ratios (in units of 10~ 7 ) in Fig. [2] 

Charm-loop effects, due to the large Wilson coefficient and the large CKM matrix element, might introduce im- 



portant effects. In a very recent work 20(, the authors have adopted QCD sum rules to investigate both factorizablc 



10 




diagrams and nonfactorizable diagrams. Their results up to the region q 2 = mju are parameterized in the following 
form, 

AC«— (q 2 ) = 1 *>° [q>q , (43) 



1 i i i )q 2 —q 
1 + r, ^^Vi 



y 9 

l + r%' ^ 

where the three results correspond to different Lorentz structures: i = 1,2,3 for terms containing V, A\ and A2 
respectively. The numerical results are quoted as follows 



AC^(f) 


- 72+ - 57 


r (l) 
' 1 


= 0.10, 


'2 " 


1.13, 




- 7fi+ - 70 


(2) 


= 0.09, 


(2) 
r 2 = 


1.12, 




- 1 11+ 114 

— L - L1 - 0.70) 


(3) 
?1 


= 0.06, 


J3) _ 
'2 — 


1.05 



(44) 

It should be pointed out that not all charm-loop effects in B — > K%l + l~ are the same as the ones in B — > K*l + l~ . 



Among various diagrams the factorizablc contributions, which can be simply incorporated into Cg given in Eq. (|A3|) . 
are the same. The nonfactorizable ones are more subtle. In particular the light-cone sum rules (LCSR) with B-mcson 
distribution amplitudes are adopted in Ref. (2o| . in which intermediate states like K* are picked up as the ground state. 
The generalization is not straightforward to the case of since in this approach states below may contribute in a 
substantial manner. However in another viewpoint, i.e. the conventional LCSR, they may be related. In our previous 
work we have shown that the light-cone distribution amplitudes of is similar with K* in the dominant region of 
the PQCD approach. If it were also the same in the conventional LCSR, one may expect that the charm-loop effects 
in the processes under scrutiny have similar behaviors with the ones in B — > K*l + l~ . Therefore as the first step to 
proceed, we will use their results to estimate the sensitivity in our following analysis and to be conservative, we use 

AC WB^K* ( _ 2) = (1 ± l)AC (i)B^' (ff 2j (45) 

in the region of lGeV 2 < q 2 < 6GeV 2 . The central values for g 2 -dependent parameters will be used for simplicity and 
in this procedure, the factorizable corrections to Cg given in Eq. (|A3|) should be set to to avoid double counting. 
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FIG. 3: Differential polarization fractions |% of B — > K^l^l - (the left panel) and B a — > f2l + l~ (the right panel). 




With the above strategy, our theoretical predictions are changed to 

f L (B^K^ + (x~) = (66.6i^)%, 
soiB^K^+fi-) = (3.49t° : ^)GcV 2 , 



f L (B s -> = (63.2±5:|)%, 

,s (i? s ^ / 2 m + aO = (3.53t°;^)GeV 2 . (46) 



The uncertainties in the zero-crossing point of FBAs are enlarged to 0.4GeV 2 . We also show the g 2 -dcpcndcnce of 
the differential polarization in Fig. [3] and the normalized forward-backward asymmetries in Fig. 21 

As a par allel way, the BR of B — > K%l + l~ can also be estimated by making use of the data of radiative B — > A'*(A" 2 )7 



decays [49 1 



B(B° -» A 2 * 7 ) = (12.4 ± 2.4) x 10~ 6 , 

S(B°^A* 7 ) = (43.3 ± 1.5) x 10" 6 . (47) 

The ratio of the above BRs R = gpftj = 0-29 ± 0.06 and the measured data of B — >• K*l + l~ shown in Tab HT1 give 
the implication 

£ OX p(5° -> Kfl + r) = (3.1 ± 0.7) x 10~ 7 , (48) 

which are remarkably consistent with our theoretical predictions within uncertainties. 

When the large energy symmetry is exploited, the seven B — > A 2 form factors can be reduced into two independent 



ones Ci and 0|- Based on these nontrivial relations, Ref. [25J has used the experimental data of B —> A 2 7 to extract 



Cl • With the assumption of a similar size for £| | , the authors also estimated the branching ratio and forward-backward 
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asymmetries of B K%l + l . Explicitly they have employed 

C ± = 0.27 ± 0.03™, 0.8U < C|| < 1.2C±, (49) 

which are comparable with our results j27| 

Cx = (0-29 ± 0.09), C|| = (0-26 ± 0.10). (50) 

As a consequence, the predicted results of BR, forward-backward asymmetries and polarizations are compatible with 
each other. 

Our results for angular coefficients, Ij = h/jj^i ar e depicted in Fig. [5] for B — > _ftT|/i + /i~ and Fig. [6] for _B S — > 
f' 2 iJ, + ljL~. Since the predictions for If, Is,Ig in the SM are typically smaller than 0.03, we shall not show them. The 
corresponding transversity asymmetries are shown in Fig. [7] and Fig. [8] respectively. One particular feature is that 
most of these results are stable against the large uncertainties from the form factors. 

For the experimental purpose, it is valuable to estimate the minimum size of the averaged value of an angular 
distribution coefficient so that it can be measured in experiment. To establish any generic asymmetry with the 
averaged value (A) of a particular decay at na level, events of the number N = n 2 /((A)) 2 should be accumulated. For 
instance on the LHCb there are 6200 events for the B — > K*l + l~ process per nominal running year JiJ. Incorporating 
all differences between K\ and K* , we may expect roughlylOOO events of B — > K%{— > Ktt)1 + 1~. Therefore if one 

wants to observe an asymmetry at no level, its averaged value should be larger than (^4) m i n = \J — 0.03n. 

Before closing this subsection, it is necessary to point out that the above estimation might be too optimistic. In 
the first few running years of LHCb, the central energy in the pp collision may not reach 14 TeV and its luminosity 
will be below 2/6 -1 . Thus in the first stage not enough data are available for a precise determination of some angular 
coefficients. Nevertheless this will not affect our analysis of branching fractions and many angular coefficients. 



B. Constraint on NP parameters and the NP effects on B — » K% A 1+ A t 



In this subsection we will first update constraints of free parameters in the above two NP models, and particularly 
we use the experimental data of b — >• sl + l~ and B —> K*l + l~ . Decay width of the inclusive process b —> sl + l~ is 
given as 



dT(b 



r/.s 



s£ + e~) I KM 2 a 2 



(i - sy< 



V cb \ 2 4tt 2 f(m c )k(m c ) 
(1 + 2s) (|C7 9 | 2 + |do| 2 ) +4(1+*) |C 7 | 2 + 12C 7 ReC. 



f(z) = 1 - 8z 2 + 8z 6 - z 8 - 24z 4 In ; 



k(z) 



2a, 
3tt 



31 

T 



a-4 2 + 



(51) 



where s = q 2 /m 2 , and T(b —> ceu e ) is used to cancel the uncertainties from the CKM matrix elements and the factor 
mj. For B — > K*l + l~, the FBAs, polarizations, and BR have been measured in different kinematic bins Q . The 



other relevant experimental data collected in Tab. [0] are from Refs. [49|, |5 



We will adopt a least-x 2 fitting method to constrain the free parameters, in which the x 2 is defined by 



the 



bt v ) 2 



(5? rr ) 2 



(52) 



where Bi denotes one generic quantity among the physical observables. The Bj he ,B^ xp and Bf r denotes the theoretical 
prediction, the cental value and 1-a error of experimental data, respectively. The total x 2 is obtained by adding the 
individual ones. It is necessary to point out that although the errors in experiment may correlate, for instance the 
measurement of B, Jl and Apg proceed at the same time in the fitting of angular distributions 0, we have not taken 
into account their correlation in our theoretical results. 

As shown in the previous section, these two NP models have the similarity that only Cc^io are modified. One 
difference lies in the coupling with the leptons, the newly introduced down-type quark in VQM will not modify the 
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FIG. 5: Angular coefficients U for B — > 



lepton sector and the coupling with leptons is SM-like; on the contrary, one new gauge boson is added in the Z' model 
and its coupling with leptons are completely unknown. 

Embedded in the VQM, the two parameters, real part and imaginary part of A s b, are found as 

ReA s6 = (0.07 ± 0.04) x 10" 3 , ImA sb = (0.09 ± 0.23) x 10~ 3 , (53) 
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FIG. 6: Similar as Fig. G3 but for B s -J- / 2 /x + M" 



from which we obtain \\ s b\ < 0.3 x 10 3 but the phase is less constrained again. The corresponding constraint on 
Wilson coefficients are 

|AC 9 | = |C 9 -C 9 5M |<0.2, |AC 10 | HC 10 -Cf M | < 2.8. (54) 
Our result of x 2 /d.o.f. in the fitting method is 49.3/(23 - 2). 
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FIG. 7: Spin amplitudes and transversity asymmetries of B — > K^fJ,^ fi 





FIG. 8: Similar with Fig.[7]but for B a -)■ f 2 fJ, + ^ 
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TABLE II: Experimental data used in the least x 2 -fitting method 



b — > civ [50J 


7 , 7-1- T — r A r\~\ 

b — > srl [49| 


T~) ^ IV * 7 -1- T — r A r\l 

B — ^ K [49| 




(10.58 ± 0.15) x 10 -2 


(3.66_q 7") x 10 


{ 1 /-\n -1-0 12 \ -in — 6 

(1.09Iq ii) x 10 




q (GeV ) 


23(10 -7 ) 


F L 


-A FB a 


[0,2] 


1.46 ±0.41 


0.29 ± 0.21 


0.47 ±0.32 


[2, 4.31 


0.86 ±0.32 


0.71 ± 0.25 


0.11 ± 0.37 


[4.3, 8.68] 


1.37 ±0.61 


0.64 ±0.25 


0.45 ±0.26 


[10.09, 12.86] 


2.24 ± 0.48 


0.17 ±0.17 


0.43 ± 0.20 


[14.18, 16] 


1.05 ±0.30 


-0.15 ±0.28 


0.70 ± 0.24 


> 16 


2.04 ±0.31 


0.12 ± 0.15 


0.66 ±0.16 


[1,6] 


1.49 ±0.47 


0.67 ±0.24 


0.26 ±0.31 



"The different convention on 81 introduces a minus sign to the forward-backward asymmetry. 



Turning to family nonuniveral Z' model in which the coupling between Z' and a lepton pair is unknown, the two 
Wilson coefficients, Cg and Cio, can be chosen as independent parameters. Assuming ACg and AC10 as real, we find 

ACg =0.88 ±0.75, AC10 =0.01 ± 0.69, (55) 

with x 2 /d.o.f. = 48.4/(23 — 2). Removal of the above assumption leads to 

ACg = -0.81 ± 1.22 + (3.05 ± 0.92)i, AC 10 = 1.00 ± 1.28 + (-3.16 ± 0.94)i (56) 

with x 2 /d.o.f . = 45.6/(23 — 4). If the /i-lepton mass is neglected, the imaginary part of C10 will not appear in the 
expressions for the differential decay widths and the polarizations. Moreover, for the forward-backward asymmetry as 
shown in Eq. (|22|) . the imaginary part of C10 contributes in the combination Rc[CgCio], thus the inclusion of Im[Cio] 
will have little effect on the x 2 ■ 

Combing the above results, we can see that the NP contributions in both cases satisfy 

|AC 9 |<3, |AC 10 |<3. (57) 

To illustrate, we choose ACg = 3 e I7r / 4 > i37r / 4 and AC10 = 3e ?7r / 4 ' l3,r / 4 as the reference points and give the plots of 
branching ratios, FBAs and the polarizations in Fig. [9] The black (solid) line denotes the SM result, while the 
dashed (blue) and thick (red) lines correspond to the modification of Cg. The dot-dashed (green) and dotted lines 
are obtained by modifying Cio- From the figure for Afb, we can see that the zero-crossing point sq can be sizably 
changed, which can be tested on the future collider or can be further constrained. 
One last process to explore is B s — > fi + , of which the branching fraction is 



B(B S -+ = TBs ^^\V; s V tb \ 2 m B J 2 ml\C 10 \ 2 I 1- ^) . (58) 



4to^ 



Using the same inputs as those in our computation of B — > K£l + l , we have 



8(B s ^fr) = 3.50 x 10- 9 ( ff; - ] ( ] %^\ . (59) 



230McV / V 4.67 



Even if Cio is enhanced by a factor of 2, the above result is still consistent with the recent measurement 51 1 

B(B S -> fi+fi') < 5.1 x 10" 8 . (60) 



VI. SUMMARY 



In this work we have explored B — > K^-^ Kn)l + l (with I = e, fj,, r) decays and a similar mode B s — > /2(1525)(— > 
K + K~)l + l~ in the standard model and two new physics scenarios: vector-like quark model and family non-universal 



17 




-0.4 1 1 I I I I I I I 1 

2 4 6 8 10 12 14 

r(QeV 2 ) 

FIG. 9: The impacts of the NP contributions on differential branching ratios (in unit of 10~ 7 ), polarization fractions and 
normalized forward-backward asymmetry of B — > K^l + l^ 



Z' model. Besides branching ratios, forward-backward asymmetries and transversity amplitudes, we have also derived 
the differential angular distributions of this decay chain. The sizable production rates lead to a promising prospective 
to observe this channel on the future experiment. 

Using the experimental data of the inclusive b — > sl + l~ and B — > K*l + l~, we have updated the constraints on 
effective Wilson coefficients and/or free parameters in these two new physics scenarios. In the VQM, we find that the 
constraint on the coupling constant is improved by a factor of 3 compared with our previous work. Their impact on 
B — > K^l + l^ is elaborated and in particular the zero-crossing point for the forward-backward asymmetry in these 
NP scenarios can sizably deviate from the SM. These results will be tested on the future hadron collider. 
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TABLE III: The values of Wilson coefficients d(mb) in the leading logarithmic approximation, with mw = 80.4GeV, /i = 
m 6 , polc [3. 



Ci 


c 2 




C4 


C B 


c 6 




1.107 


-0.248 


-0.011 


-0.026 


-0.007 


-0.031 


-0.313 4.344 -4.669 



Appendix A: Effective Hamiltonian 



The effective Hamiltonian governing b — > sl + l is given by 

10 



s = —^v tb v*^2c i ( f ,)o i (n), 

v 2 i=l 



(Al) 



where Vtf, = 0.999176 and Vj s = —0.03972 |50j are the CKM matrix elements and Cj(/x) are Wilson coefficients for the 
effective operators Oi. In this paper, we will adopt the Wilson coefficients up to the leading logarithmic accuracy |48|, 
and their values in SM are listed in Tab. IIII1 Since the NP scenarios considered in the present paper would not 
introduce any new operator, the SM operators will form a complete basis for our analysis 



Oi 
3 



= {s a Ca)v-A(cpbp)v-A) 2 = (s a Cp) V -A{C[]b a ) V -A, 



V-A, 



{s a b a )v-A'^2{qpqi3)v-A, Oi = (s a b/3)v-A22(qpq a ] 

5 = (s Q 6 Q )y_A^(q/3(?^)y+A, 6 = (Sab^y-A^iqpqaW+A, 



o 9 



~2T 



-±sa^(l- l5 )bF^, 



(h,l)(s^(l - 75)6), O10 = ^7,750(^(1 ~ 75)6) 



(A2) 



The left-handed and right-handed operators are (gi92)y-A(Q3'Z4)v'±A = ( | ?i7 A '(l-~75)92)(<?37^(l±75)'3'4)- "ib = 4.8GcV 
and m s = 0.095GeV are b and s quark masses in the MS scheme and a cm = 1/137 is fine structure constant. The 
double Cabibbo suppressed terms, proportional to V u bV* s , have been neglected. 

At the one-loop level accuracy, the matrix element of b — > sl + l~ transition receives loop contributions from 0\ — O^. 
Since the factorizable loop terms 1521 can be incorporated into the Wilson coefficients C7 and C9, it is convenient to 
define combinations Cf ff and C| ff |52J 



Cf 
Cf{q 2 ) 



C7 — C5/3 — c 6 , 

din) + h(m c , s)C - s)(4C 3 + 4C 4 + 3C 5 + C 6 ) 
-h(Q, s)(C 3 + 3C\) + ^(3C 3 + C 4 + 3C 5 + C 6 ), 



(A3) 



with s = q 2 /ml, Co = C\ + 3C*2 + 3C3 + C4 + 3Cs + 6*6, and rh c = m c jmh- The auxiliary functions used above are 



/»(*,$) = --m^--lnz+- + -z--(2 + :r)|l-c C | 1 / 2 



In 



l-x-l 

2arctan- 



i7r for ir 



< 1 



for x = —r- > 1 



/i(0,s) 



8 . m b 
■— in 

9 /Lt 



9 



■ Ins 



27 



(A4) 



In the following, we shall also drop the superscripts for C| and C~f s for convenience. 

On the hadron level resonant states, such as vector charmonia generated from the b — > ccs, may annihilate into a 
lepton pair. Therefore they will also contribute in a long distance manner But these contributions can be 

subtracted with a kinematic cutoff in experiment. Moreover our following analysis of differential distributions will be 
mainly dedicated to the region of lGcV 2 < q 2 < 6GeV 2 , also excluding contributions from the charmonia. 
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Appendix B: Helicity amplitudes 



Within a graphic picture B — ¥ K 2 {-^ Kir)l + l~ proceeds via three steps: B meson first decays into an onshcll 
strange meson plus a pair of lcptons; the K 2 meson propagates followed by its strong decay into Kir. To evaluate the 
decay width of multibody decays, we shall adopt the helicity amplitude which mainly uses 

g, v = -J2^W<W + ^- (Bl) 

A q 

e is the polarization vector with the momentum q and A denotes the three kinds of polarizations. The last term can 
be formally identified as a timelike polarization e A1 (t) = -%=, and thus the metric tensor <7„„ can be then understood 

as summations of the four polarizations. For the purpose of illustration we will first evaluate the decay amplitude of 
B —> K 2 l + l~ ■ In the SM, the lepton pair in the final state is produced via an off-shell photon, a Z boson or some 
hadronic vector mesons. These states may have different couplings but they share many commonalities: the Lorcntz 
structure for the vertex of the lepton pair is either V — A or V + A or some combination of them. Therefore the decay 
amplitudes of B — > K 2 l + l~ can be rewritten as 

A(B^K* 2 l+r) = &(L)Hp{V)+L*(RyHp{R\ (B2) 

in which C ^(L) , C ^(R) are the lepton pair spinor products: 

C^L) = h^l - 75 )l, C^R) = h,(l + 75K (B3) 

while % incorporates the remaining B — > K 2 part. In the case of masslcss leptons, left-handed and right-handed 
sectors decouple, which will greatly simplify the analysis. The identity in Eq. (|B1[) results in a factorization of decay 
amplitudes 

A(B -> K*l+r ) = C^{L)U v {L)g^ + C^R)U v {R)g^ 

= -^ClxHlx-^CrxHrx, (B4) 

A A 

where is the momentum of the lepton pair and Cl\ = ^{L)e^(X) and £r\ = £ A1 (i?)e M (A) denote Lorcntz invariant 
amplitudes for the lepton part. It is also similar for the Lorentz invariant hadronic amplitudes: Rlx = R^i^e^X) 
and Rr\ = H M (i?)e* (A). The timelike polarization gives vanishing contributions in the case of m; = for I = e,/j,: 
using equation of motion, this term is proportional to the lepton mass. 

An advantage of the helicity amplitudes is that both hadronic amplitudes and leptonic amplitudes are Lorentz 
invariant. Such a good property allows to choose different frames in the evaluation. For instance leptonic amplitudes 
are evaluated in the lepton pair central mass frame, while hadronic B decay amplitudes are directly obtained in the 
B rest frame. Since K\ and K* have several important similarities, B — > K 2 l + l" differential decay widths can be 
simply obtained from the ones of B K*l + l~ in a comparative manner. 

• Longitudinal and transverse B decay amplitudes are obtained by multiplying the factor ^ n f^ m an d 

respectively. The function A is the magnitude of the K% momentum in B meson rest frame: 

A = \(m 2 B ,m 2 K ,,q 2 ) = 2mB\pK*\, and A(a 2 ,6 2 ,c 2 ) = (a 2 — b 2 — c 2 ) 2 — 46 2 c 2 . This replacement is an out- 
put of the fact that the polarization vector e is replaced by ct in the form factor definitions. Explicitly, these 
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hadronic amplitudes are 



H 



LO 



N 



\f%m B m K * 2m K *J7f- 



(Cg - C w )[(m% - m 2 K , - q 2 )(m, B + m K *)Ai 
A 



m B + m K * 



+2m b (C 7L - C 7R )[{m B + 3m 2 K * - q 2 )T 2 - — — j 3 



m „ — rn 



Kr, 



H 



L± 



= N 



m B m K * 



(Cg - C w )[(m B + m K :)Ai =p 



2m b (C 7L + C 7R ) ^ ±% /^ T ^ + 2m b (C 7L - C 7fl ) ^ 



(m| - mi*)T 2 



(Cg — Cio) — p=^4o, 



\/8m B mT 
H L i\c 10 ^-C 10 



(B5) 



withiV = -^2 r F tb V t *. 

• In the propagation of the intermediate strange meson, the width effect of K% could be more important since 
T K * ~ lOOMeV > T K , ~ 50MeV Nevertheless, since the ^ width is only larger than that of K* by a 
factor of 2, the narrow- width approximation, which has been well used in the case of K* , might also work for 
K 2 . In this sense, there is no difference except that the B(K* — > Kir) is replaced by B(K 2 — ^ Kw). 

• Incorporation of the K% — > Kir decay gives the complete results for differential decay distribution of B — > K%(—t 
Ktt)1 + 1~. Angular distributions of K 2 and K* strong decays are described by spherical harmonic functions: 



1 ' 



for K* and Y 2 (6, (ft) for K 2 . In particular we find the relations 
~3 

~3 

~7T 



cos(9 K ) = C{K*) 
■sm(9 K ) = S{K*)- 



16tt 
15" 



(3 cos 2 9 K - 1) = C(K%), 



32tt 



sm{29 K ) = S{K^). 



(B6) 



Our formulas for branching fractions and forward-backward asymmetries can be shown compatible with the ones 
in Ref. [25j through the following relations 



Alo = N K *a L m B 



1 



A L ± 
A, 



2m K * V q 2 

2\N K *B T 

-{A-£), 



(1 - mi. - f)T + XQ + 0-- rh 2 - f)B -AC), 



2m B 
2\N K ,p T 



2m B 
N K *a L VA 



(B7) 

(B8) 
(B9) 
(BIO) 



m K * yJq A 

where the coefficients A,B,8,J-,Q,T-L are defined in Eqs. (49,50,53, 54) in Ref. but the coefficient C in Eq. (51) 
contains a typo and should be read as 



C = 



1 



1 ™ m 2 K . 



(1 - m K ^)cf{s)A^ (a) + 2m 6 cf (if 2 (s) + 



ff / rriK'2 



1 — m'L* ^* 

— ^rf 2 ( s )) 



The dimcnsionless constants are given as A = A/m B , to_jst* = mx* /m B , rh b = m b /m B and q 2 = q 2 /m 2 B . oil 
and p T = l/x/2. 



(Bll) 
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